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Abstract 

The coupling between the intramolecular vibrational modes and the doped 
conduction electrons in M^Cqq is studied by a calculation of the electronic 
contributions to the phonon self energies. The calculations are carried out 
for an orientationally ordered reference solid with symmetry Fm3m and for 
a model with quenched orientational disorder on the fullerene sites. In both 
cases, the dispersion and symmetry of the renormalized modes is governed by 
the electronic contributions. The current current correlation functions and 
frequency dependent conductivity through the midinfrared are calculated for 
both models. In the disordered structures, the renormalized modes derived 
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from even parity intramolecular phonons are resonant with the dipole excited 
single particle spectrum, and modulate the predicted midinfrared conductiv- 
ity. The spectra for this coupled system are calculated for several recently 
proposed microscopic models for the electron phonon coupling, and a com- 
parison is made with recent experimental data which demonstrate this effect. 

PACS Numbers : 78.20.-e, 74.70 Ya, 74.20.-z 

I. INTRODUCTION 

Soon after the discovery of solid phases of C 60 M it was demonstrated that these solids 
can be doped from an insulating parent phase through a series of interesting conducting 
and nonconducting phases by the introduction of alkali species in the interstitial volume 
0. There has been particular interest in the properties of the family of M 3 Cqq (where 
M is an alkali metal) which exhibit superconducting ground states. The interplay of the 
superconductivity and the unique structural properties exhibited by these molecular metals 
is a key issue for these materials. 

Among the various proposed theories for the superconductivity of these phases is a model 
which posits that an effective attractive interaction between the conduction electrons is me- 
diated by exchange of intramolecular phonons in the doped solid ^M- Since this model was 
originally proposed, it has been appreciated that the theory can be tested quantitatively by 
examining the renormalization of the intramolecular vibrational modes which derive from 
their coupling to the doped conduction electrons. Indeed, the key Raman active vibrational 
modes are experimentally found to shift and broaden in the conducting state confirm- 
ing the underlying idea that coupling between the intramolecular modes and the conduction 
charge is important for the physics of the doped phases. 

This theoretical picture is complicated in an essential way by the orientational "merohe- 
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dral" disorder || of the alkali doped samples which have been synthesized. Experimentally, 
the structure of most of the M 3 Cqq samples appears to be that of an orientational glass. In 
these solids, each of the Ceo molecules, on the lattice sites of an FCC Bravais lattice, adopts 
one of tho high symmetry orientations with the two fold molecular axes oriented along 
the (001) crystal directions. There are two inequivalent ways in which each molecule can 
achieve this orientation, and they are related by a 7r/2 rotation about any (001) direction. 
This choice of orientation varies, essentially randomly, from site to site in the solid. 

This structural disorder has important consequences for virtually all the electronic prop- 
erties of the doped solid, including the renormalization of the intramolecular modes men- 
tioned above. Nevertheless, the effects of the disorder are awkward to treat properly in the 
theory, and consequently most theoretical work has treated the effects of orientational disor- 
der only very approximately. Indeed, a considerable body of work has investigated the "zero 
band width" limit of this problem 0,§]||, in which the intermolecular hopping amplitudes 
are set to zero and progress can be made analytically. This approximation is based on the 
fact that the intermolecular electronic hopping amplitudes are much smaller than the in- 
tramolecular amplitudes which allows a useful separation of energy scales. Thus many solid 
state properties do reflect the underlying molecular structure of the solid. For example, the 
conduction band relevant to our discussion below is derived primarily from the antibonding 
ti u molecular state ||10|| , with very little mixing between different intramolecular orbitals in 
the solid phase. Similarly, each of the vibrational bands in the solid may be identified with 
a fundamental vibrational excitation of the isolated Cgo- This intramolecular symmetry also 
strongly constrains the possible couplings between the vibrational degrees of freedom and 
the conduction electrons : electronic density fluctuations on a given fullerene site may be 
coupled only to the fully symmetric a g or quadrupolar h g vibrational excitations on the same 
site §. 

Of course the zero band width limit does not address any phenomenon which derives 
from the intermolecular properties of the doped solid. For example, the two point density 
density response function for the doped carriers depends explicitly on the three dimensional 
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dynamics of the doped carriers in these molecular solids. This in turn depends sensitively 
on the degree and type of microscopic structural disorder in the doped phases, so that 
spectroscopy on the low energy electronic excitations can be used to probe the structural 



disorder in the doped solids [TT|,[T2|. In this theory, scattering from static orientational 
fluctuations in the doped solid leads to a residual mean free path for the conduction charge 
of order 20 A |]13 |, and to a residual resistivity of order 3000 S/cm JL^]. This strong scattering 



is believed to dominate the observed low frequency conductivity which is still lower, of order 
1500 S/cm ||l4| , [l5l at low temperatures. 



The phonon renormalizations mentioned previously are obtained from a similar family of 
two point functions, and are also quite sensitive to the intersite orientational correlations of 
the solid. For a reference orient at ionally ordered solid, the high symmetry of the structure 
leads to selection rules which are applicable to the solid phase and which further restrict 
the coupling of the vibrational and electronic degrees of freedom. These selection rules can 
be violated in the disordered solid, depending on the type and degree of the orientational 
disorder. This is very important for understanding the renormalization of the intramolecular 
vibrations in both the ordered and disordered phases, and the response functions derived 
from these modes flTB|. In fact, the spectrososcopic signatures of this effect are now beginning 



to be reported in experiments for the doped fullerides where several of the even parity "on 
ball" vibrations can be observed as modulations in the dipole active long wavelength and 
low frequency electronic excitation spectra. This is strictly forbidden for the orientationally 
ordered structure, and its appearance in experiments provides a powerful direct probe of 
the scattering by oreintational fluctuations in the doped phases. This effect is reminscent 
of similar phenomena which have been previously encountered in the charge transfer salts 
T7HI9[ , and in the doped conducting polymers [^Df , although the structural origins of the 



symmetry breaking here are somewhat unique to the fullerides. 

In this paper, we present a study of the electron phonon coupling in the doped solids in 
which we explicitly treat the three dimensional dynamics of the doped conduction charge. 
This theory is carried out first for a reference orientationally ordered structure, and then for 



a model with orientational randomness. Below we calculate the electronic contributions to 
the phonon self energies in both models, so that the effects on orientational disorder on the 
phonon self energies can be clearly identified. Finally, we compute the phonon contributions 
to the long wavelength current current response function, and study the infrared conductivity 
for this coupled electron phonon system. A brief comparison with experiments is given at 
the end of the paper. There we find that the mid-infrared conductivity for the coupled 
system depends sensitively on the disorder and the distribution of electron phonon coupling 
strengths, and may therefore be used to critically compare several popular models for the 
electron phonon coupling constants for these materials. 



II. DISORDER AND ELECTRONIC STATES IN M S C, 



3^60 



In this section we review the tight binding description of conduction electronic states in 



the doped fullerides. Since this construction has been presented in detail elsewhere |pTT] , |2T 
our presentation will be brief. 

The conduction band of the alkali doped fullerides is derived from an antibonding molec- 
ular orbital of the isolated Cqq with t lu symmetry. This is the LUMO of the neutral molecule; 
it is three fold degenerate, and transforms like a vector under spatial rotations. The band 
structure of the doped and undoped orientationally ordered phases has been studied by first 
principles methods p2|- p^1 . For both cases it is found that the conduction bands can be 
described reasonably well by a model restricted to hopping of a three component field = 
(cj/i,/^ = x,y, z) between nearest neighbor fullerene sites: 

H ei = Y,Yl T ^ ) c lc> ( 1 ) 

where ej„ creates an electron at site i in the /xth orbital. The matrix elements in T depend on 
the orientation of the bond between sites i and j, on the orbital polarizations of the two sites 
connected in the hop, and on the molecular orientations of the terminal sites. In accordance 
with experiment || we will restrict our attention to the situation where the terminal sites 
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can adopt either of the two high symmetry settings with the two fold axes aligned along the 
(100) crystal directions, which we label the "A" and "B" orientations. The choice of setting 
can be indexed by an Ising spin variable, s, defined on each site and which takes the value 
Si = 1(-1) for the A(B) orientation. The intermolecular hopping amplitude then depends on 
the spin variables on neighboring sites in the manner, 



T^) = T(°\n J )+T^(T lJ )s l + T^(r^ 



(2) 



+T^(r ij )s l s J 



So that the hopping amplitude can be developed from the four fundamental hopping matrices 
T^. For a hop along a reference bond flj = (a/2)(i +j)=f the matrix T take the form 



T®(f )=t 



T^\r )=t 



ABO 
B A 
c 

X -Y 
Y —X 




T^(r ) = t 



T&(f )=t 



X Y 
-Y -X 


DEO 
EDO 
OOF 



(3) 



where the three basis orbitals are assigned to the x, y, and z orbital polarizations, respec- 
tively. The hopping amplitudes for any other bond direction can be obtained by a rotation 
of the T's listed in equation 3. The values of the parameters A, B, C, D, E and F , deduced 
from the work of Gelfand and Lu, are listed in table 1. We should note that our convention 
for assigning the orbital polarizations is fixed in the crystal frame and therefore does not 
change when the molecular orientation is rotated. This is in fact different from the conven- 
tion originally adopted by Gelfand and Lu, and the reader is directed to the Appendix A of 
Ref. [^TJ for a discussion of this point. 

Below we will investigate the effects of orientational disorder in the supercell approxima- 
tion in which we construct a periodically repeated unit cell containing a large (typically 32) 
molecules with the orientations randomly assigned to the A and B settings. The electronic 
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Hamiltonian for this model is developed from the fundamental hopping matrices listed in 
equation |3]. The artificial periodic boundary conditions allow us to analyze this problem in 
the Bloch orbital basis : 



C k,i)i ^/jy e Ic I,ifi (4) 

in a solid with N supercells, indexed by the supercell translation vectors Ri, and where i and 
\i index the basis site and orbital polarization respectively. In this basis the Hamiltonian for 
a supercell containing M sites is a 3M x 3M matrix H(k) , for which the eigenvectors C nk 
are expanded in the Bloch orbital basis: 

Cnk = ^A n ^{k)c\ i ^ (5) 



Several authors [pTl^ [T3[| have applied this model to study the conduction electron states 
in orient at ionally ordered and disordered structural models. The fluctuations in intersite 
terms in H e i due to the disorder in this problem are on the order of the Fermi energy, and 
so one can correctly anticipate that the effects of disorder will be quite significant. Explicit 
calculations demonstrate that the full conduction bandwidth in the disordered phase is only 
slighly smaller than the bandwidth obtained for the orientationally ordered structure and in 
addition the density of states at the Fermi energy is not significantly different from the value 
for the ordered structure. However, the electrons are strongly scattered by the orientational 
fluctuations, and the calculated conduction band spectrum for the disordered alloy is smooth 
and relatively structureless. A closer inspection of the underlying conduction electron states 
in the disordered orientational alloy reveals that electron propagation near the Fermi level 
can be well described in terms of an effective Bloch wavevector k, and a mean free path 



of order 20 A fI5| . There is no evidence in the theory for localization of the one electron 
states in this model except in a narrow range of energy near the band extrema which is 
irrelvant for electronic processes near the Fermi energy. This is in accord with experimental 
measurements which show a metallic conductivity at low temperature in the normal phase 



15| , that is the measured or extrapolated conductivity at low temperature comfortably 



exceeds the minimum metallic conductivities expected at these electron densities. Thus 
theory and experiment provide a consistent overall picture that the electronic structure near 
the Fermi energy for the disordered alloy is that of a "dirty" metal with propagating states 
at the Fermi energy and a relatively short residual mean free path. 



III. INTRAMOLECULAR PHONONS 

In this section we present the model for the "bare" unrenormalized phonons of the doped 
solid. The 176 modes of a single fullerene molecule can be indexed using the ten irreducible 



representations of the icosahedral group Ih [ 25| , |26 |. However, the conduction states are 



derived from the t\ u molecular states, and therefore the conduction charge in this manifold 
may couple only to the lattice modes of a g and h g symmetry. There are two such fully 
symmetric modes, they are the molecular breathing mode and pentagonal pinching mode 
of the icosahedral molecule. There are eight families of h g modes, and these describe the 
various possible quadrupolar deformations of the Cqq molecule. These intramolecular modes 
span an energy range from 200 cm~ l to 1600 cm~ l . 

In our model, the bare phonon Hamiltonian will be represented by the Einstein spectrum: 

H ph = Y,Y, Uam(bl m (i)b am (i) + -) (6) 

i am 

where &„ m («) creates the m-th mode of symmetry a on the i-th site. We shall use a con- 
vention where a = 0(with m = 1,2) denotes the a g modes and a = l,5(with m = 1,8) 
label the h g quintuplets. The degeneracy of each h g quintuplet implies that for a = 1,5, 
^am = u m . These even parity modes are experimentally accessible by Raman scattering, 
and measurements of the Raman frequencies for M x Cqq with x = 0, 3 and 6 have now been 
reported from several groups . The striking result of the experiments is that the Raman 
lines in M^Cqo are shifted and broadened relative to the undoped solid. On the other hand in 
MqCqq, which is an insulator, although the frequencies are further shifted, the Raman lines 
are again sharp. It is suggested that the source of the shift and broadening in the M 3 C 60 
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phase is primarily the coupling of the intramolecular modes of equation 6 to the conduction 
charge. For the MqCqq structure which is a band insulator, the charge gap is of order 0.5 
eV, so that the phonon self energy is again purely real, leading to a shift of the frequency 
but no damping. 

There are several important contributions to the frequency shifts and linewidths of the 
Raman active vibrations of the doped system. Most notably, the intramolecular interatomic 
potential is intrinsically anharmonic. Thus addition of charge, even to an isolated molecule, 
will lead to a relaxation of the bond lengths, and to a shift of the bare frequency. Here 
"bare" implies that the phenomenon is not a solid state, but rather a molecular effect. Note 
that this effect is not accompanied by electronic induced broadening which is a solid state 
effect. The frequency renormalization in the M 6 C 60 insulating phase is due to just such an 
"intramolecular" effect. 

The situation in the interesting M^Cqq phases is complicated due to the confluence of 
these "intramolecular" and "intermolecular" solid state processes. The added conduction 
charge induces an intramolecular relaxation of the bond lengths and frequency shifts of 
the type observed for the M 6 C 6 o phases pP|J2E|f27[ ] . There are, in addition, self energy 



contributions due to the "one loop" density fluctuations of the added conduction electrons. 
(These are clearly abent for the terminal phase with x=0 which has no conduction charge, 
and for the x=6 phase which has a saturated band). These terms are nevetheless the 
most interesting since they communicate information about the solid state structure, and 
particularly about the intersite orientational correlations, into our effective Hamiltonian for 
the lattice degrees of freedom. These terms are clearly essential for what follows as they 
induce a mixing, frequency shift and damping of the coupled intramolecular modes. 

For simplicity we shall assume that the two types of phonon renormalizations, inter- 
molecular (due to the one loop corrections from the conduction charge) and intramolecular 
(due to the anharmonic terms in the bare Hamiltonian and to higer order terms in the elec- 
tron phonon coupling) are simply additive. Thus the bare Hamiltonian of equation |^ will be 
assumed to already contain the intramolecular effects, and we will explicitly calculate the 
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intermolecular contributions absent from that model. To be precise, Faulhaber et. al. |28 
have reported theoretical results for the lattice vibrations of the isolated Cfj" trianion, which 
properly contain the intramolecular contributions only, and we adopt their results for the 
frequencies to an appearing in equation |[ For convenience, these frequencies are reproduced 
in Table 2. We then proceed to explicitly evaluate the important "solid state" one loop 
electronic contributions to our effective Hamiltonian in the doped phases. 



IV. ELECTRON PHONON INTERACTION 

The lattice modes of equation || are coupled linearly to the conduction electrons via an 
on site deformation potential 

#e^ = EEE r >^A« (7) 

i A V-v 

where 

Qx(i) = b x (i) + b{(i) (8) 

is the displacement corresponding to the A — th mode on site i. Here V* is a matrix acting 
in the electronic t\ u subspace. The symmetry of this matrix is fixed solely by the symmetry 
of the coupled mode A while the strength of the interaction depends on the particular mode 
of interest. If we write A = am for the m-th mode of symmetry type a we have that 

where a = denotes a fully symmetric a g mode and a = (1,..,5) denote the five basis 
functions in the h g representation of Ih- Anticipating our passage to a crystal with cubic 
symmetry we will break this manifold into the invariant e g and t<i g subspaces and write: 



a g 



73 
73 / 
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A 1 



75 







A 2 



^0 0^ 

75 
V° 75 



A 3 







V 75 





f i 





] 




V2 





°] 




e a : A 4 = 





1 





A 5 = 





1 





(10) 
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v 
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We have defined these matrices so that they are orthonormal with trA Q A^ = 5 a! 3. As 
expected, the form of the A matrices requires that the symmetric modes couple linearly to 
the total on ball conduction charge density and the h g modes are linearly coupled to the 
electric quadrupole density fluctuations. 

The (? Qm 's appearing in equation 9 are the deformation potentials are obtained by study- 
ing the linear variation of the intramolecular Hamiltonian with respect to the normal coor- 
dinates of the vibrational mode of interest. 

n 



9ctm A 



[IV 



/ i dH mo i | . 
" (/i| dQ~~ W) 



(11) 



where H mo i is the intramolecular electronic Hamiltonian and M is the mass of a carbon atom. 
A microscopic evaluation of these coupling coefficients requires a model for the dependence 
of H mo i on intramolecular displacements, as well as a microscopic assignment of the various 
intramolecular phonon frequencies and eigenvectors. There have been several evaluations 
of the g a nS reported in the literature ||f|,|28| which treat H mo i and the bare vibrational 
modes at various levels of sophistication. The coefficients g an reported in these studies are 
displayed in table 3. Note that for h g modes (a = 1, 5), g am = g m . We also note that the 
data tabulated here differ from those reported in several of the original papers owing to 
different normalization conditions on the A matrices adopted by the various authors. These 
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various models may be classified according to two groups: those that localize the strongest 
interactions to the two highest frequency h g modes and those that distribute the coupling 
more uniformly over the phonon spectrum. The two highest frequency h g modes contain 
double bond stretching character, and therefore might be expected to couple strongly to the 
added conduction charge which occupies an antibonding state. The coupling to the lower 
frequency modes is a more subtle matter, deriving from the perturbation to H mo i due to 
bond angle fluctuations. 

The expression for H e i- V h is valid for an isolated molecule. Information about the solid 
state properties is then contained in the electronic amplitudes (cj^c^) appearing in this ex- 
pression. Using the notation in Section II for a periodically repeated supercell, the linearized 
interaction reads 

Hei- Ph = E E E r^cuMi) (12) 

Transforming to the basis of Bloch eigenstates C\ k given by the expression, 

C\ k = Y. A nUkY k - Rl c\^ (13) 

Iifi 

we obtain, 

Hd-ph= E ^nn'(Kk')CUk)C n ^k')Q Xi (k-k') (14) 

k,n,k',n' Xi 

where, 

QM=b Xi (-q) + bl i (Q) (15) 

and , 

f^(*, = E ^Ki,{k')A nliv {k) (16) 

is the electron-phonon coupling evaluated in the Bloch basis. The phonons relevant to our 
discussion of the infrared absorption spectra are those with q=0 and hence we will set k = 
k' in the above expressions and abbreviate the diagonal components T^ n ,(k, k) by T^ n ,(k) 
and Qm(0) by Q Xi . 
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V. PHONON SELF ENERGIES 



In this section we calculate the one loop corrections to the dynamical matrix for the 
phonons in the x=3 phase. The calculations are carried out first for a reference orienta- 
tionally ordered model, and then for a model containing quenched disorder. We focus on 
the selection rules governing the electron phonon coupling for the ordered case, and the 
breakdown of these rules in the orientationally disordered phase. We also compute the 
renormalized phonon Green's function which includes the effects of the one loop corrections 
and which will be used to calculate the current current correlation function for the coupled 
electron phonon system in Section VI (b). 

We will study the phonon Green's function 

D(\i, AY; u) = -i dte- l ^(T[Q Xt , Q[ H ,]) (17) 
Jo 

As in the previous section, A indexes the a particular symmetrized on ball displacement 
field, and i indexes a site in the unit cell. We specialize the discussion to the q=0 part of the 
Green's function, and therefore enforce periodic boundary conditions on the displacement 
field Qxi+Qxi- In the absence of the electron phonon interaction the unrenormalized Green's 
function is diagonal in both the site and phonon "band" indices A. 

D (Xi, Xi'; u) = D (X, Lu)5 X \'bn> (18) 

where 

D (X,u)= 2 f Ux , (19) 

This Green's function is assumed to contain all the relevant intramolecular renormalizations, 
which are then parameterized by the bare frequencies uj\ and the linewidths e\. For the 
former we use the theoretical data of Ref. |p8|| , while for the latter we use the experimentally 



observed linewidths for undoped C 60 reported in |7|]. These parameters are listed in Table 
3. 
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The linear electron phonon interaction in equation |7| is represented by the elementary 
vertex shown in Figure 1(a). In the Bloch orbital basis, this matrix element is diagonal 
in the site index but in general not in the orbital indices fi. This scattering amplitude 
provides a second order contribution to the phonon self energy, where the relevant term 
describes a process whereby an electron hole pair, created on a common site % by the vertex 
shown in Figure 1(a) propagate and ultimately recombine on a site %' in the solid. The 
vibrational degrees of freedom on different sites are therefore coupled indirectly through the 
density fluctuations of the conduction electrons. This process, represented diagramatically 
in Figure 1(b) , is evaluated to give the self energy contribution 



U(\i,Xi ;uj) = EfcE„^ 



a E n u — E n ih — to — i5 



nk J-'n'k 

(f(E nk ) - f{E n , k )) (20) 

For this molecular solid in which the "bare" phonons are well described by Einstein oscillators 
with negligible direct intersite coupling the indirect coupling in equation |20] determines the 
dispersion of the vibrational bands of the doped phase. Similarly, all the information about 
the orientational structure of the solid phase enters our theory through the spatial variation 
of the one electronic states which appear in the polarization contribution in equation ED. 



The separable form of the deformation potential in equation ^ allows us to write the 
self-energy term in equation in the factorized form: 

n(Az, AY; u) = g am IL(ai, pi'; u)gp m > (21) 

where a labels the phonon symmetry, and n is a phonon band index. 

For the orient at ionally ordered phase we can drop the site index since we have only one 
molecule per unit cell. In addition since the icosahedral symmetry is still maintained, the 
reduced polarization term II is diagonal in the symmetry labels a, so that 

U{a,i3) = 8 aP U{a) (22) 

Furthermore, since we have chosen the coupling matrices (A) to be irreducible under cubic 
symmetry we have the additional simplification, 
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n(o) ± n(i) = n(2) = n(3) ± n(4) = n(5) (23) 



Ofcourse, the full II in equation |2T] still has off-diagonal terms in the phonon band indices. 
Thus, in the orientationally ordered solid, the electron phonon interaction induces coupling 
only between phonons of the same symmetry in different bands. 

However, in the orientationally disordered phases, off diagonal terms can and do appear 
even in II and they induce a significant mixing not only between phonons with different 
symmetries but also between phonons on different sites. In addition, there is also a significant 
variation in the the diagonal components of II from site to site as we see below. 

We note that in equation ^0] we have omitted the long wavelength contributions from 
the "intra-band" terms with n = n'. In fact, for coupling to the (optical) Einstein modes 
relevant to our model, in which uo — > constant as q — > 0, these contributions vanish. This is 
in striking contrast to the situation for coupling to acoustic modes, in which an additional 



intraband contribution also exists. Allen [|31]] has derived an approximate expression for 



Im(II) for general wavevector q which is of the form, 

Im(n(g)) = u £ \g(k, k + q)\ 2 6(E k )6(E k+q ) (24) 
k 

The q = limit of this expression is indeed non-zero and results in a frequency shift (Acu n ) 
and linewidth(7 n ), 

2 2 

Au n = N(0)^; ln = nco n N(0) 2 ^ (25) 

where d n is the degeneracy of the nth mode and iV(0) is the density of states at the fermi 
level. Although this result has been previously applied to the doped fullerides [0,0], it 



is strictly applicable only in the limit that the Debye frequency for the relevant acoustic 
phonons is much smaller than the fermi energy, Ep, for the conduction charge. These 
conditions are violated for the fullerides where we have optical phonons with frequencies 
comparable to the fermi energy. In our treatment we will therefore consider only the self 
energy terms appearing in equation 22. 
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Once the self energy contributions are calculated using equation |20| the renormalized 
phonon Green's function is obtained by solving the Dyson equation 

D = D Q + D Q UD (26) 

For the orient at ionally ordered solid, the factorization in equation (|2TJ) allows us to reduce 
the matrix equation (pq ) to a scalar problem in each invariant subspace. For the disordered 
solid in which mixing between the various symmetry branches and sites occurs one has to 
solve the full matrix equation. However, the mixing between the h g and a g manifolds is 
expected to be small and we solve the equation in each subspace. Due to the factorisation in 
(pl|) the largest matrix which needs to be inverted for the h g phonons is only 5M x 5M (where 
M is the number of sites). 

For the orientationally ordered phase, the reduced polarization terms for the fully sym- 
metric space a = vanish identically. Thus a long wavelength excitation in the fully sym- 
metric subspace provides a constant shift to the electronic potential, but cannot scatter 
electrons between conduction states. The remaining nonzero contributions in the irreducible 
e g and t 2g subspaces are displayed as functions of the phonon frequency in Figure 2. The 
imaginary part of the self energy (shown in the top panel) introduces a damping of the 
molecular modes, the real part (lower panel) describes the frequency shift due to the screen- 
ing. We see that the overall shape of Im II is similar in the two subspaces, peaking near 
2000 cmT 1 . We find that the t 2g vibrations are coupled slighly more strongly, particularly 
at higher frequency. The most important consequence of this is that the real part of the self 
energy at low frequency ( below 1000 cmT 1 ) is enhanced in the t 2g channel relative to the 
e g channel. This implies a "crystal field " splitting of the £=2 intramolecular modes with 
the t 2g frequency systematically suppressed relative to the e g frequency. Solving equation 
(25) we find that this splitting is of order 30cm" 1 for the second h g mode( h g (2) in Table. 
2) whose renormalised frequency is 372 cm -1 . Finally, we note that the higher frequency 
Raman active modes overlap the region in which Im II is large. This leads to a relatively 
large damping of the higher frequency Raman active vibrations. However, as we see below 
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this situation is complicated somewhat by the orientational disorder of the host. 

For the orientationally disordered phase the symmetry selection rules, (equations ([22| 
and |23p no longer pertain and there are a large number of independent components of the 
reduced self-energy, EL In Figure 3 we display a representative sample corresponding to 
the diagonal components of EI corresponding to the same symmetry branch but different 
sites(i.e. U(ia = 5,i/3 = 5) for various values of site i). It can be seen that there is 
significant variation in II on different sites. This spread is expected to lead to an additional 
inhomogenous broadening of the phonon linewidths. This shape of these spectra can be 
understood simply. If we model the conduction band by a rectangular density of states and 
assume a contact interaction between the electrons and phonons which is independent of 
frequency and wavevector, the predicted Im EI (a;) is obtained from a simple convolution of 
the filled and empty sectors of the one electronic spectrum. This yields a density fluctuation 
spectrum which is symmetric with a triangular shape, as indeed we find from our numerical 
computations. 

The analogous self-energy terms in the a g subspace are shown in Figure 4. The shape is 
again "triangular" and the overall scale is similar to that found for the h g subspace. Note 
however, that the self energy contributions are quite different owing to the variation of the 
amplitudes g a . It is nonetheless significant that these contributions are nonzero for the 
disordered phase. This reflects the fact that the long wavelength phonons in the folded 
Brillouin zone of the supercell are coupled to intracell charge fluctuations, and can thus 
be screened and damped by their coupling conduction electrons. The simplest example of 
this nonzero coupling in the symmetric subspace is the branch which is backfolded from q 
= (^,0,0) in the extended zone, where n is the supercell period. This couples to the x- 
polarized density fluctuation of the conduction charge in the supercell. 

We note that there are two intramolecular modes of a g symmetry: the breathing mode 
near 200 cm -1 and the pentagonal pinch near 1500 cm" 1 . These occur respectively well 
below, and near, the peak of Im EI(u;) shown in Figure 4. Of course the phonon self energy 
depends on both the magnitude of fl(ou) and the coupling strength g. The coupling constants 
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for the a g modes were computed in Ref . j| ( see Table 3) where it was found that the coupling 
to the breathing mode was very small while the coupling to the pentagonal pinch mode is 
quite substantial. The combination of this large coupling strength with our calculated II for 
the a g symmetry channel would predict a substantial broadening of the higher frequency a g 
mode. On the other hand experimentally ]7J the increase in the linewidth of this phonon is 
small which seems to indicate that the actual coupling should also be small. 

VI. INFRARED CONDUCTIVITY 

In this section we use the renormalized phonon Green's function developed in the previous 
section to calculate the frequency dependent conductivity of the coupled electron phonon 
system. We begin by considering the purely electronic contributions to the current current 
correlation function for the ordered and disordered phases. We then turn to the indirect 
contributions to the conductivity mediated by coupling through the dressed phonons. 

A. Electronic Contributions 

The real part of the frequency dependent conductivity is obtained within linear response 
theory as the imaginary part of the fourier transform of the current current correlation 
function. Since the wavelength of the probe is long on the scale of atomic dimensions it 
suffices to study the q = limit of these response functions. In our tight binding formalism, 
the paramagnetic term in the current operator is given by 

3 = ie rtiTp, ( 27 ) 

ij V 

where T is the intermolecular hopping matrix along the T^-th bond. Thus the paramagnetic 
current-current correlation function is calculated from the Kubo formula 

Xab(u)=i dte-^([j a (t)J b (0)]) (28) 

J 
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where the angular brackets denote a thermodynamic average, and a, b are Cartesian indices. 
The frequency dependent conductivity, which is the linearized current in response to an 
external electric field is therefore 

CTab{U) = (29) 

For the noninteracting electronic Hamiltonian introduced in Section II, this correlation func- 
tion can be readily calculated from the single particle eigenstates to yield 

(<7ei)a&M= jfEk^nn'3nn'( k )jn'n( k ) ( 30 ) 

8{E nk -E n , k -u){ f{ Y k -E E * k) ) 

where f(E) is the Fermi-Dirac distribution, and jn n >(k) is the matrix element of the a-th 
component of the current operator in (26) between single particle eigenstates n and n'. 
In this treatment which neglects the diamagnetic contributions to j, the long wavelength 
low frequency Drude response needs to be evaluated with some care. One may proceed by 
studying the intraband contributions to j at finite wavevector q, and extract the q — > limit: 

df k 



«r») = £JLWJL(£0^%) = ds( u ) (3i) 

kn 0tj nk 

so that the weight of the Drude pole is obtained as a Fermi surface average 

D = Y,fnn(k)j b nn(k)S(Enk) (32) 
kn 

Rather than compute the Fermi surface average required in (30) our strategy will be to 
develop the sum rule for the total conductivity 

- / doo(a e ) ab (uj)) = e 2 (K ab ) (33) 

7T JO 

where, in the tight binding formalism (K ab ) is the ground state expectation value of the 
operator: 

Kab = [Ra,jb] = J2z\,( T ij)^ T ij)bT,,u(ij)cl u C ju (34) 

ij V 
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Equation |34] is the generalized long wavelength f-sum rule p2] within our tight binding 



model. For the conducting phase, the interband paramagnetic part of the response function 



in |3T] does not exhaust this sum rule, so that the residual weight provides the spectral weight 
in the Drude pole. For the free electron gas, the interband terms vanish identically, and then 
full spectral weight is then located in the pole at zero frequency. Our strategy is therefore 
to evaluate the sum rule, and the interband response function directly, and to extract the 
Drude piece by a comparision of their respective integrated spectral weights. 

This method can be applied to both the orientationally ordered and disordered phases. 
The results for the frequency dependent conductivity are displayed in Figure 5, and the 
corresponding partitioning of the interband and "intraband" spectral weights are given in 
Table 4. For the orientationally ordered phase we observe a broad midinfrared continuum, 
extending from zero frequency to 4000 cm _1 peaking near 1600 cm' 1 . This peak can be 
roughly identified with a peak in the joint density of states do to the two primary peaks 
in the single particle density of states for the ordered phase. Approximately 40 % of the 
total spectra weight is located in the Drude pole. We also remark that the interband 
response function extends to zero frequency, where it vanishes linearly with u. This is 
due to an interesting symmetry in the band structure for the orientationally ordered doped 
phase: the Fermi energy crosses a doubly degenerate band along the (111) directions in 
reciprocal space. This degeneracy is lifted as on shifts away from this symmetry direction, 
yielding a continuum of low frequency interband excitations which extends exactly to zero 
frequency. We remark that this degeneracy has some interesting consequences for the theory 
of superconductivity on such a Fermi surface which we have discussed in a previous paper 



Our results for a supercell containing quenched orientational disorder ( also see |12[) are 
shown in the lower panel in Figure 5. The supercell for this calculation contains 32 sites. 
Here we see that the interband peak is broadened and shifts to lower frequency where it 
overlaps the overdamped remnant of the Drude pole. The composite spectrum thus exhibits 
an overdamped midinfrared response function with a distinctly non-Drude lineshape. We 
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remark that due to the artificial periodic boundary conditions in our supercell model, this 
spectrum retains a very small but nonzero spectral weight in a true Drude pole at zero 
frequency. This weight is nevertheless very small, and accounts for only 2 % of the total 
spectral weight in this disordered model. This nonzero residual weight is nevertheless an 
artifact of the periodic boundary conditions, and will vanish as the size of the disordered 
cell is increased. 

The results of Figure 5 demonstrate that merohedral disorder has a very strong effect 
on the dynamics of the conduction electrons, a point which has been emphasized in several 
previous studies [pTT|— jT3| . The predicted low frequency low temperature conductivity is of 
order 3000 S/cm. This is rather small, but significantly larger than the experimental low 



temperature value of 1500 S/cm ||14| . We believe that the discrepancy arises, at least in 
part, due to the effects of direct electron electron interactions which have been omitted from 
our treatment. 



B. Lattice Contributions 

Here we discuss corrections to the direct current current correlation function due to the 
coupling of the conduction charge to the intramolecular phonons. We are not concerned with 
the odd parity intramolecular vibrations which may be directly coupled to a long wavelength 
dipole active probe. Rather we will consider the even parity intramolecular modes which 
are indirectly coupled through the response conduction electrons. The appearance of these 
modes in the infrared is a very interesting phenomenon which, as we noted previously, can 
not occur for the orientationally ordered structure, nor indeed any for structure which pre- 
serves the inversion symmetry. Thus the appearance of these features in experimentally 
measured spectra is a fundamental signature of the microscopic breaking of inversion sym- 
metry by orientational disorder. As we see below, these spectra are also quite sensitive 
to the details of the electron phonon coupling model, so that experimental measurements 
of this effect might be used as a critical test of various microscopic models. The relevant 
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microscopic process which indirectly couples the even parity intramolecular modes to light 
is illustrated diagramatically in Figure 6. As with the construction of the "direct" current 
current correlation function, the external field is directly coupled to density fluctuations of 
the conduction charge. These may in turn be coupled to the intramolecular modes through 
the vertex of Figure 1. The composite polarization term is one which thus correlates the 
current operator with the vibrational deformation potential and will be referred to as the 
"phonon charge." Thus the dressed phonon Green's function is indirectly coupled to the 
external field through this effective phonon charge, contributing to the full current current 
correlation function of the coupled electron phonon system. The microscopic theory of this 
process has been studied for low dimensional CDW systems jlTlJl8| and for the conducting 
polymers p0| . In addition, it has been recently employed to explain the dopant induced 
infrared activity of the t\ u vibrations of the Cqq ions in the A x Cqq compounds p4| . As we 



noted above, the process represented in Figure 6 must vanish identically in an inversion sym- 
metric structure. To explicity demonstrate this, it is conventient to examine the nonlocal 
conductivity for the coupled system in the position representation: 

Xph(w) = J2 v ( Xi ~> u + iS)D{\i, X'j; u + i8)V(X'j; u - iS) (35) 

AA' 

where D is the renormalized phonon Green's function , V(Xi) is the electronically induced 
phonon charge for the A -th phonon on the i-th site. V(Xi) is given by 

V(Xt; = [ dE^\j ilvl G m {E)T Xi G w {E + u)\ (36) 

where G is a 3x3 matrix one particle Green's function in the space of orbital polarizations. 
In equation (35) the subscripts i are site indices and the trace is taken over the orbital 
indices. The operator T Al describes the intramolecular scattering of a conduction electron 
due to its linear coupling to the A-th mode on site i: 

^ iiu = r c^Cjyj (37) 

Since the isolated fullerence molecule is inversion symmetric, it is clear that T Xl is invariant 
under the inversion operation, P. 
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p—l^atp ^£ai (38) 

Similarly, the current operator j appearing in equation should transform as a vector, which 
is odd under P. 

P-'jP = -3 (39) 

We note that on the lattice, the operator j involves the "off diagonal" elements of the elec- 
tronic hopping operator T and thus equation (38) pertains only in the presence of inversion 
symmetry. Finally, for the inversion symmetric structure the one particle Green's function, 
and the two point density response function derived from it is even under parity. Thus the 
phonon charge, which involves a trace over a product of these three terms is odd under 
parity, and therefore vanishes in the q — > limit. For the disordered system, equation (36) 
is still true. However, the effective current operator in (26) is no longer odd under P, and 
furthermore, the two point density response function does not commute with P. In this case 
the effective phonon charge for an intramolecular "even parity" phonon is nonvanishing. 

We now turn to an evaluation of the of indirect contribution to the conductivity given 
in equation (34). Our disordered system is modelled by a supercell, with periodic boundary 
conditions on the molecular orientations. In such a model 

Xph (uj)= ]T V{Xi;u + iS)D(Xi,X'i , ;u + iS)V(Xi , ;u-iS) (40) 

Xi,\'i' 

where now i and %' are site indices, and the phonon charge is explicitly: 

V(Xl] U + l5) = f Efc Enn' p p ^ 

(f(E nk ) - f(E n , k )) (41) 

Note that in this model, crystal momentum is conserved modulo a reciprocal lattice vector 
of the supercell, and equations (|40|) and ( f4"T|) therefore represent the appropriate q=0 limits 
of these response functions. 

The effective phonon charge in equation ([Hp is evaluated directly by numerical inte- 
gration on a mesh of 12 3 wavevectors in the folded Brillouin zone reciprocal to a supercell 
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containing 32 molecular sites. The calculated contributions to Re a p h due to coupling to the 
h g modes are displayed in Figure 7, where we have collected results for three different models 
for the electron phonon interaction. These are obtained for each of the coupling amplitudes 
g a listed in Table 3. The contributions of equations fl4"0|) and ([II]) to the conductivity are 
additive and therefore each of these spectra should be superposed on the direct contribution 
to the conductivity shown in Figure 5(b). We see that for all the models, the coupling to the 
intramolecular "quadrupolar" modes produces a substantial modulation of the otherwise 
smooth featureless midinfrared conductivity. These modulations are correlated with the 
frequencies of the renormalized h g phonons, although it is found that the lineshapes are typ- 
ically relatively broad and quite complex. The overall scale of these modulations is of order 
100 S/cm for all the models studied. Interestingly we observe that the detailed structure in 
the modulated midinfrared conductivity is different for the various models. As we remarked 
earlier, the electron phonon models can be broadly grouped into two categories. In the first 
of these, the electron phonon coupling is dominated by the coupling to the two highest lying 
h g branches of the spectrum. The model due to Varma et. al. is a well known example 
from this group. These two highest lying modes are dominated by double bond stretching 
character on the fullerene, which are expected to couple strongly to the " antibonding" t\ u 
electronic state. In the second group of models, the coupling is distributed more uniformly 
among both the high frequency and lower frequency h g vibrations. The models derived by 
Schluter et. al. and by Faulhaber et. al. belong to this second category. These models 
require substantial coupling to the bond angle fluctuations and the conduction electrons. 
We find that in these latter models, but not in the former, the indirect contributions to the 
conductivtiy shows two well defined peaks in the range between 300 and 500 cm -1 . In this 
region, the imaginary part of the phonon self energy is relatively small, and therefore the 
midfrequency intramolecular modes, if coupled to the conduction charge, show up as sharp 
features in the spectra. By contrast, for the former class of electron phonon models, with 
negligible coupling to the midfrequency modes, the sharp low frequency features are absent 
and the modulation spectra are "rugged" over a rather wide energy range, extending to the 
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highest frequency modes nearl700 cm-1. 

The process illustrated in Figure 6, also leads to an indirect contribution to the conductiv- 
ity from the high frequency, fully symmetric a g phonon branch. This resulting contribution 
to the conductivity is displayed in Figure 8 calculated using the theoretically calculated 
value from Ref. || Here the antisymmetric Fano type lineshape which is the signature of 
the resonance of this local mode with the pair continuum is clearly evident. The scale of 
this modulation is of order 100 S/cm, and its width is due to the relatively large coupling 
constant g a for the second a g mode listed in Table 3. Experimentally, this mode is observed 
as a much narrower resonance, so that the microscopic prediction for the coupling strength 
for this mode appears to be overestimated. 

Modulations of the midinfrared conductivity in K^Cqq have recently been reported in 
experiments by Kuzmany et. al. |35| , |36| In Figure (9) we reproduce from this work, the real 
part of the full midinfrared conductivity, and the indirect phonon contribution, which is ob- 
tained by subtracting a smooth electronic background from the experimental measurements. 
We note that the experimental data in this Figure is taken at T = 300K and the scale of 
the conductivity over this range of frequency is known to depend on temperature. With 
this caveat, the modulations of Re a due to coupling through the phonons is clearly evident 
in the experimental spectra. The two very sharp features near 600 cm and 1400 cm are 
odd parity modes, which are "directly" dipole active and thus not described by the indirect 
phonon charge model derived above. The remainder of the modulations in the conductivity 
are on a scale of typically 25-50 S/cm, and an additional sharp feature due to coupling to 
the high frequency symmetric a g mode is clearly evident. These modulations provide mi- 
croscopic spectroscopic evidence for the local breakdown of inversion symmetry due to the 
orientational disorder in these samples. We also note that the experiment does not show 
any sharp phonon structure in the midfrequency range. Thus we believe that the type I 
electron phonon coupling models provide a more faithful representation of the experimental 
situation in these fullerides than the type II models. 
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VII. SUMMARY 



In this paper we have focused on the interference between the intramolecular vibrational 
modes and the low lying single particle excitations measured in the midinfrared conductivity. 
Recent experimental data do provide some supporting evidence for this coupling in the doped 
phases. However there are three aspects of the measured spectra which remain relatively 
poorly described by this theory, and to which we should call attention. First, the overall 
scale of the conductivity is incorrectly predicted in the theory. In fact the calculations 
overestimate the low frequency limit of Re u(uj) by nearly a factor of two. This is a serious 
difficulty in that the predicted value of a(uo — > 0) is not rescaled by a simple rescaling of 
the conduction bandwidth. We believe that the failure to predict the correct scale of the 
measured conductivity is related to a deeper problem, and is likely associated with the neglect 
of direct electron electron interactions in the model. Second, careful experiments at very 
low frequency show a rather sharp minimum in the frequency dependent conductivity near 
lu ~ 100cm -1 . This is thought to separate a low frequency Drude like peak, from a rather 
broad midinfrared continuum. This sharp low frequency minimum is absent from the theory 
presented above, even in the presence of coupling to the phonons. It is possible that this 
originates in some strong couping to intermolecular vibrations neglected in our treatment. 
Finally, our discussion does not address the very interesting temperature dependence of the 
resistivity, which increases roughly linearly with temperature over a very broad temperature 
range, extending to 700 K. At the highest temperatures studied, one may infer from these 
data that the effective mean free path is smaller than the radius of the C 60 , a state of 
affairs which appears difficult to reconcile with the theoretical picture given above. Thus, 
while we believe that the phonon structure in the midinfared conductivity can be adequately 
explained by the theory presented here, there are a number of issues to be resolved in the 
experimental low frequency spectra. 
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Table 1 Hopping Amplitudes for the electronic Hamiltonian(Equation [1] obtained from 

0.014 meV for doped 



Ref. |PH| . These coefficients are multiplied by a scale factor,^ 
fullerides. 



A 


B 


C 


D 


E 


F 


X 


Y 


0.01 


0.38 


-2.29 


2.09 


-2.36 


0.38 


-0.63 


-0.49 



Table 2. Frequencies and line- widths for the bare phonons( in cm 1 ). The frequencies 
are as calculated in Ref. |28| and the linewidths are chosen to agree with experiment ( 0). 



mode 


Frequency 


Linewidth 


h 9 (l) 


244 


4.2 


h a (2) 


412 


5.5 


h g (3) 


560 


7.5 




750 


9.0 


h a (5) 


1073 


7.0 


h g (6) 


1202 


7.6 


h 9 (7) 


1395 


7.5 


h g (8) 


1508 


9.5 


a g (l) 


507 


2.5 


a g (2) 


1470 


1.5 



Table 3. Coupling const ant s,g an , as defined in Equation [Ll|, for the h g and a g modes(in 
meV) 
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mode 



h g (2) 
h g (3) 

h g (4) 
h g (5) 
h g (6) 
hg(7) 
h g (S) 



%(2) 



Ref. j 

9.4 

7.1 

8.9 

0.0 
29.0 

0.0 

77.8 

46.0 



Ref. | 
11.6 
13.5 
27.6 
20.5 
0.0 
0.0 
26.5 
44.0 



0.0 
52.0 



Ref. |28j 

6.0 
21.1 

7.6 
28.7 
13.4 
18.0 
47.8 
37.8 



Table 4 Spectral weights in units of te 2 /ha (where a is the lattice constant) correspond- 
ing to the interband and total electronic conductivities. 



Model 


Interband weight 


Total weight 


Ordered Crystal 
Disordered Crystal 


15.2 
26.6 


25.0 
27.1 
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FIGURE CAPTIONS 



FIG.l. (a) The electro n-phonon intereaction vertex in the Bloch orbital basis. (b) The 
one loop contribution to the phonon self energy,!!. The electron and phonon propagators 
are represented by a solid and wavy line respectively. 

FIG. 2. Imaginary(upper panel) and Real(lower panel) parts of the reduced phonon self- 
energy ( 11(a)) in the orient at ionally ordered solid corresponding to the e g ( a = 1,2) and 
t2g{ ot = 3,4,5) symmetry channels. 

FIG. 3. Imaginary (upper panel) and Real(lower panel) parts of the diagonal compo- 
nents of the reduced phonon self energy (Jl(ia,ia)) for h 9 phonons in the disordered solid 
corresponding to a = 5 and a representative sample of sites i. 

FIG. 4. Imaginary (upper panel) and Real(lower panel) parts of the diagonal compo- 
nents of the reduced phonon self energy (Il(ia,ia)) for a 9 phonons in the disordered solid 
corresponding to a = and a representative sample of sites i. 

FIG. 5. Real part of the pure electronic conductivity for the ordered(upper panel) and 
disordered (lower panel) solids 

FIG. 6. Diagrammatic representation of the process which generates the phonon contri- 
bution to the conductivity. The solid line denotes the electron propagator, the wavy line is 
the renormalised phonon propagator and the dotted line the photon propagator. 

FIG. 7. Results for the contribution of the h 9 modes to the conductivity using various 
models for the coupling constants (g an ) listed in Table 3. 

FIG. 8. Results for the contribution of the high frequency a g mode to the conductivity 
using the coupling calculated in Ref. || (see table 3). 

FIG. 9. (a) The experimental mid-infrared conductivity |35| with a fit representing the 
smooth electronic background (b) The contribution due to the phonons obtained by sub- 
tracting the electronic background. 
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